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ABSTRACT 

An experimental rotating wickless heat pipe apparatus was 
designed and machined. The apparatus includes a rotating heat pipe 
assembly, test stand, spray cooling assembly, safety shielding, 
aga instrumentation. 

A revised condensing limit for the operation of the rotating 
heat pipe was obtained by modifying Ballback's Nusselt film con- 
densation theory to include the effects of a thermal resistance in 
the condenser wall and in the condenser outside surface cooling 
mechanism. Approximate results, obtained for half-cone angles of 
1, 2, and 3 degrees, show that less heat can be removed than originally 
predicted by Ballback[7 ], and that the outside heat transfer co- 
efficient can significantly alter the condensing limit. 

An improved Nusselt theory was developed which applies for 
all half-cone angles, and which includes the effects of the thermal 
fepistances in the condenser wall and in the condenser outside surface 
cooling mechanism. This formulation led to a second-order non-linear 
differential equation for the film thickness which was numerically 
integrated using a free-overfall boundary condition at the condenser 
exit. Results obtained for a half-cone angle of O degrees are sub- 
stantially less than the results obtained from the approximate solution 


for half-cone angles of 1, 2, and 3 degrees. 
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I. INTRODUCTION 


A. CONVENTIONAL HEAT PIPE 

A heat pipe is a self-contained device capable of transporting 
large quantities of heat at nearly isothermal conditions. Con- 
ventionally it consists of three main parts: a container, a working 
fluid, and a wicking structure capable of promoting capillary action. 

When heat is added to one end of the container, some of the 
working fluid evaporates. This end is the evaporator. The pressure 
differential in the vapor chamber due to the increase in vapor at 
the evaporator causes the vapor to flow to the opposite end of the 
container, transporting heat as latent heat of vaporization. Here 
the vapor condenses on the cooler surface and releases the heat to 
be removed from Che structure. This end of the container becomes 
the cont Coed The condensate is pumped back to the evaporator 
section by capillary action, thus completing the cycle. The heat pipe 
is unique in that it operates with no moving parts, maintains nearly 
isothermal conditions, and removes the dependency on gravity through 
the use of capillary pumping. 

Limits on the heat transfer papab ladies of conventional heat 
pipes are imposed by a number of fluid dynamic mechanisms; including a 
wick resistance limit, a sonic vapor velocity limit, a vapor velocity 
entrainment limit, and also : wick boiling limit. These limits have 
been considered in a number of studies [1, 2, 3, 4 ]* with the con- 
clusion that in general a conventional heat pipe using an ordinary 
working fluid is limited by the rate that the condensate can be returned 
to the evaporator due to wick resistance. 


* Numbers in brackets indicate references listed in Bibliograpry- 
ou 








B THE ROTATING, WICKLESS HEAT PIPE 

The restriction imposed by the inability of the wicking 
structure to provide sufficient pumping force to return condensate 
to the evaporator can be alleviated by removing the wick and 
utilizing centrifugal acceleration as a pumping source [2] Lane 
pipe is allowed to rotate about its longitudinal axis and the wall 
is given a slight taper, a rotational force field can easily be 
produced that not only overcomes the viscous losses in the system, 
but also the losses due to a hindering gravitational Ii1eld] sine 


operation of a rotating heat pipe is shown schematically in FIGURE 1. 


Heat In Heat Out 


bes te stbotahenles Lt td] 





Liquid Annulus Condensate Return 

BOILER SECTION CONDENSER SECTION 
FIGURE 1. Schematic Drawing of Typical Rotating Heat Pipe 

The applications of a rotating heat pipe appear to be wide and 
varied. Wherever there is rotation and a need to transfer heat, the 
use of such a system is feasible. The cooling of motor rotors, turbine 
rotor blades, aircraft braking systems, and air-conditioning units are 
just a few examples of the uses to which the rotating heat pipe concept 


can be applied. 
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The capillarity of the wicking structure need no longer be 
considered as a limit on heat pipe operation since it has been 
replaced as a pumping force by an artificial gravity head. Ballback 
[7] has studied the limitations imposed on the rotating system by 
various fluid dynamic mechanisms in a preliminary theoretical analysis 
of the operation of a ee heat pipe. By using existing theories 
and experimental correlations, he was able to estimate the limits on 
rotating heat pipe operation imposed by the onset of critical nucleate 
boiling, condensate entrainment, and choked vapor flow. An estimate 
of a condensing limit was obtained by performing a Nusselt-type analysis 
for condensation on a truncated cone, since existing theory and 
correlations were not applicable. This analysis, however, was restricted 
by certain simplifying assumptions that were made to obtain a closed- 
form analytical result. Ballback's results indicated that the operation 
of a meee ine heat pipe is limited by the amount of heat that can be 


removed from the condenser section. 


C. THESIS OBJECTIVES 
The objectives of this study were to: a.) design and construct 
an experimental rotating heat pipe assembly; b.) revise and supplement 


the initial theoretical study of rotating heat pipe operation performed 


by Ballback. 
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II. EXPERIMENTAL PROGRAM 


A. BASIC DESIGN CONCEPTS 

For Ballback to generate numerical results from his theoretical 
model, a fixed geometry for the heat pipe had to be decided upon. 
He chose to model the evaporator section alone the lines of the rotating 
boiler used at the Lewis Research Center in Cleveland, Ohio, for the 
study of boiling heat-transfer coefficients operating under high gravity 
levels [6]. With this geometry fixed he allowed the size of his 
condenser section (length and base radius) to be governed by the amount 
of taper in the section. The taper was defined by the half-cone angle, 
which was a parameter in the theoretical study. A strong influence in 
the overall design was the desire to be able to visually observe the 
mechanisms that were taking place within the heat pipe during operation. 
The onset of peice nucleate boiling and the entrainment of the 
returning condensate were two principal phenomena which were to be 


observed. 


B. DESCRIPTION OF EQUIPMENT 

The components of the wickless heat pipe are divided into four 
sub-groups to simplify the explanation of their design and construction. 
These are: 1.) evaporator 2.) condenser, 3.) auxiliary equipment, and 
4.) instrumentation. FIGURE 2 is a schematic representation of all the 
components of the test apparatus. A cross-sectional drawing of the 
assembled heat pipe is pictured in FIGURE 3 and FIGURE 4 is a photograph 


of the machined pieces of the rotating heat pipe prior to assembly. 
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1.) The Evaporator 
The evaporator was modeled to approximate the rotating 

boiler tested at Lewis Research Center [6] . It was machined of 
stainless-steel because of this metal's corrosion-resistant properties. 
It has an inside diameter of 3.125-inches, an exit diameter of Paee 5 - 
inches, and an overall length of 5.905-inches (see FIGURE 3). This 
length includes 2/3-inch for the insertion of two pyrex-glass view- 
ports, 4-inches of heating surface, and an intended $-inch adiabatic 
section. The inside diameter of the evaporator is stepped 4-inch from 
the exit diameter to allow for the formation of a liquid annulus. The 
inner viewport rests on a teflon-coated metallic O-ring seal and is 
separated from the outer viewport by a 3/32-inch thick compressed-fiber 
gasket. The viewports are held in place by a stainless-steel end cap 
which is secured in place by twelve +-inch stud-bolts. These bolts 
also support two bronze collector rings. The rings are insulated from 
their supports by 1/16-inch teflon tubing and washers. Electrical power 
is passed by a brush assembly consisting of eight spring-loaded graphite 
brushes through the collector rings to an ll-gage Chromel-A heater wire 
in a 3/16-inch inconel sheath. Power to the heater is furnished by a 
DC motor-generator capable of delivering 150 amperes at 220 volts. The 
heater was helically wound and silver-soldered in place to insure good 
thermal contact. After the heater was coated with a thin layer of 
Sauereisen cement, asbestos insulation was packed around the outside 
to reduce the radial heat losses. In addition two radial 1/16-inch 
grooves were machined on the outside surface to within 1/16-inch of 
the inside boiler wall at either end of the heating element to lessen 


axial heat losses. 
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With this construction and by assuming a maximum operating 
temperature for the heater element, it is possible to estimate the 
maximum heat flux into the boiler. Likewise, it is possible to 
estimate the heat losses by approximating the outside convection 
coefficient from an analysis of air flow over a stationary cylinder 
[11 }. These calculations are shown in Appendix A. The results 
indicate that a maximum radial heat flux of approximately 470,000 
Btu/nr-ft* can be obtained with an associated radial heat loss of 
approximately 1300 Benner (with 3/4-inch asbestos insulation). 

As noted in Appendix A, the heat flux is dependent upon the 
thermal resistances that lie in the radial heat transfer path. The 
dominant resistance is that of the stainless-steel evaporator wall 
because of its low thermal conductivity. It is therefore advantageous 
to keep the thickness of the wall as small as possible. However, as 
mentioned in the section on instrumentation, thermocouples are to be 
imbedded in the boiler wall to monitor the temperature to obtain a 
temperature gradient and associated heat flux. From an uncertainty 
analysis (Appendix B), it is found that the farther apart the thermocouples 
are radially mounted to one another, the smaller the uncertainty in 
determining an accurate heat flux. So there is a trade-off on the 
selection of the wall thickness. In this case, a thickness of ?-inch 
was used. 

A photograph of the machined evaporator section is shown in 


FIGURE 5. 
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2.) The Condenser 

The condenser section was machined of stainless-steel in 
the shape of a truncated cone with a 3 degree half-cone angle. 

The base of the condenser mates to the evaporator exit and has an 
inside diameter of 2.5-inches (see FIGURE 3). A step is created 

at the evaporator-—condenser intersection where the condensate flows 
from the condenser to the recessed inside diameter of the evaporator. 
The condenser tapers from an inside diameter of 2.5-inches to one of 
1.46-inches at the truncated end over a length of 10-inches. From here 
the condenser extends an additional 1 1/2-inches to house a cylindrical 
end plug. Of the overall length, 9-inches (along the axis) is con- 
sidered to be the cooling surface. The first inch of the condenser 
surface near the evaporator is assumed to be adiabatic. The wall 
thickness of the condenser is 1/16-inch to keep the thermal resistance 
as small as possible. 

The condenser section has two flanges. The large front flange 
mates to the bearing-scat flange on the evaporator. They are held 
together by the previously-mentioned 1/4-inch stud-bolts. A flange at 
the small end mates to a similar flange on the end plug, through which 
there are eight 1/8-inch stud-bolts for assembly purposes. The end 
plug has an extention on its open end that has four l-inch square 
surfaces upon which are mounted four phenolic junction boards. 

A photograph of the condenser and end plug (with shaft in place) 
is shown in FIGURE 6. 

3.) Auxiliary Equipment 

The auxiliary equipment are sub-divided into four groups 
for ease of description. These are: a.) drive assembly, b.) test 
stand, c.) spray cooling assembly, and d.) safety shields. 
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a. Drive Assembly 

Te drive assembly consists of the end-plug, shaft, 
pulley, drive motor, and bearing. The end plug has a 3/4-inch 
threaded shaft way machined into it (See FIGURE 3). Between the 
shaftway and end of the end plug is a .27 inch hole to allow for the 
escape of non-condensable gases during ptlling (to be described later) 
and for the insertion of a pressure transducer into the vapor space 
during operation. The shaft is a stainless-steel cylinder which is 
hollowed to allow the passage of thermocouple leads from the junction 
poards to a slip-ring unit. ‘The shaft has 3/4-inch of its outside 
diameter threaded for screwing into the shaftway, and is threaded for 
3/8-inch along its inside diameter on both ends. The condenser end 
holds a #-inch arm that extends through the hole in the end plug into 
the vapor space, and into which iS screwed 2 small, semi-conductor pressure 
transducer, i. 9 teflon-covered thread for sealing. The opposite 
end of the shaft is threaded a accomodate a coupling to the slip-ring 
unit that rotates with the pipe. 

One-half of a square keyway is machined on the shart and 
one-half on the shaftway such that they coincide when the shaft is 
screwed securely into place in the end plug. A stainless-steel square 
key is slid into the keyway. The key transmits the rotational torque 
from the shaft to the heat pipe. The torque is applied to the shaft 
via a 2.65-inch pulley and v-belt by a 2-horsepower, 4-phase variable 
speed motor. The motor has a maximum speed of 4500 rpm and a minimum 
speed of 450 rpm. The speed is controlled by a push-button electric 
remote control unit. The motor is equipped with a magnetic disc brake 


for rapid shut-down. 
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A double row, angular contact bearing is mounted on a 
seat machined on the shaft. This bearing and the large (5.96-inch 
inside diameter), floating, single bearing (with bakelite compound 
cage to withstand high rpm) mounted on the bearing seat on the evap- 
orator flange allow the pipe to easily rotate in any orientation. 

be Test Stand 

The test stand was designed to sturdily support the 
heat pipe assembly mounted on a 15x44-inch steel bed-plate and the 
motor on its bed-plate. The system operates in orientations ranging 
from 0 to 90 degrees from the horizontal. Both bed-plates Sy eueneen 
to a 4-inch iron pipe. The pipe is placed in three- 2-inch thick 
clamps which are supported 2-feet off the ground by +-inch steel plate 
As the pipe is rotated to any orientation, the motor rotates with it 
and the assembly is then clamped securely in place. Vibrations are 
kept to a minimum by vibration-mounting the motor. FIGURE 7 is a 


drawing of a side and front view of the test stand configuration. 


i 
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Heat Pipe Motor 
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FIGURE 7 Schematic Drawing of Test Stand Configuration 
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Ce Spray Cooling Assembly 
The cooling of the condenser section is accomplished 

by spray cooling. The asSembly consists of a cylinder rolled from 
1/8-inch thick stainless-steel with an outside diameter of 13-inches 
that encloses the condenser section. The cylinder is cut in half 
longitudinally, The bottom half is secured to the bed-plate on which 
the heat pipe is mounted and the other half is bolted to its mate 
once the heat pipe is in Place and ready for operation. Each half cylinder 
holds four spray nozzles at the same axial Pesmrron, ~ Te spray nozzles 
are easily interchanged so that mists of droplets ranging in size 
from 300 to 600 microns can be obtained. To keep the spray from 
interfering with the operation of the heater or either of the bearings, 
1/8-inch thick shields are welded to the ends of the half cylinders. 
When the cylinders are mated these Shields come within 1/16-inch of 
the bottom ae Sides of ridges machined on both the front and rear 
condenser flanges (see FIGURE 3). 

The coolant flows from copper feed lines to a steel mixing tube 
mounted on top of the spray cylinder. Here the average temperature 
of the coolant is measured by a quartz thermometer. - From the mixing 
tube the coolant flows through plastie tubing to each of the nozzles. 
Once sprayed it passes through a drain in the cylinder and bed-plate 
to a second mixing tube which is also monitored by a quartz thermometer. 
From here the coolant is sent to an external drain. 

d.) Safety Shields 
The entire heat pipe is surrounded by 1/8-inch 


stainless-steel Shielding to insure safe Operation. The shielding consists 
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of three components: the evaporator housing, condenser nousing, 
and shaft housing. 

The evaporator housing is an inverted-U shaped section 
that encloses the evaporator. The front end is closed and contains 
a plexi-glass window to observe internal pipe operation through 
the pipe viewports. The back end is open nee mates with the spray 
cooling cylinder end shields. The siding has a curved plexiglass 
window for the visual observance of exterior pipe operation and a 
door located near the bed-plate for the adjustment of the brush 
assembly. The inverted-y section slides onto set screws on the bed- 
plate for easy mounting. 

The condenser spray cooling cylinder serves not only as 
the spray cooling assembly but also as the condenser housing safety 
eageld. 

The shaft housing is also shaped like an inverted U and 
encloses the shaft, pulley, and rear bearing assemblies. Both ends 
are left open with one sliding against the spray cooling cylinder 
end shields. A slot is left in the siding for the operation of the 
v-belt. This section also slides onto set screws on the bed-plate. 

4. Instrumentation 

The instrumentation is grouped into the following sub- 
divisions for ease of explanation: a.) temperature measurements, 
b.) pressure measurements, c.) coolant flow rate measurements, 

d.) rotational speed measurements, and e.) electrical power measurements. 


A schematic drawing of the instrumentation setup is shown in FIGURE 2. 
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a.) Temperature Measurements 

Temperature readings are obtained from thermocouples 
mounted in the evaporator wall to determine a temperature profile 
and associated heat flux through the wall. Four grounded copper- 
constantan thermocouples with a 1/16-inch diameter stainless-steel 
Sheath are utilized. They are placed Pee eee diameter 
wells drilled 33-inches deep into the stainless-steel evaporator. 
To insure that no minute air gaps are present between the well walls 
and the sheaths, a coating of Silver Goop is placed on the thermo- 
couple prior to insertion. Two of the thermocouples are mounted 
opposite each other at the same radial position to check on the 
uniformity of the temperature field. This pair is mounted at the 
middle of the }-‘nch wall thickness. The other two thermocouples 
are also mounted opposite one another, one 1/16-inch radially from 
the outside evaporator suface, one 1/16-inch radially from the inside 
Sutrace. 

The leads from the thermocouples are brought through both 
evaporator and forward condenser flanges, down the side of the condenser, 
and through the rear condenser flange to phenolic junction boards. 

From the boards the leads go inside the shaft and down to the slip-ring 
unit which is attached to the main assembly by a coupling mounted at 
tne end of the shaft. 

A thermocouple is suspended in the liquid annulus of the evaporator 
to measure the bulk liquid temperature during operation. A sixth 
thermocouple is suspended in the vapor space at the evaporator exit to 
record the saturation temperature of the vapor. The leads to the 


thermocouples in the liquid and vapor pass through slots in the evaporator 
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wall (see photograph in FIGURE 4) after passing through the condenser 
forward flange. 

The temperature of the coolant is measured before and after spraying 
by quartz thermometers mounted in mixing tubes located above and below 
the condenser spray cooling cylinder. 

b.) Pressure Measurements 

The saturation pressure of the vapor is measured by a miniature 
semi-conductor pressure transducer mounted in the vapor space. The 
transducer is temperature compensated and is threaded on a ¢-inch 
Giameter arm that extends from the shaft through the end plug and 
into the vapor space. The probe's leads pass throuch the arm, up 
through the shaft to the junction boards, back into the shaft. and 
down to the slip-rings. 

The pressure of the coolant prior to spraying is indicated by 
a pressure tap on the coolant feed line. 

ec.) Coolant Flow Rete Measurements 

The coolant flow rate is Ace by passing the coolant 
through a flow rotameter prior to spraying. 

ad.) Rotational Speed Measurements 

The rotational speed of the heat pipe is found by two different 
methods. First a Hewlett-Packard optical tachometer is employed. Six 
equally spaced spots mounted circumferentially on the shaft cause the 
tachometer to emit an electrical impulse each time one passes through 


its optical field. These impulses are fed to a Donner frequency 
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counter and counted over a ten second interval. The read-out is in 
revolutions per minute. The second method of determining the rpm 
is by the use of an electronic strobe light. 
e.) Electrical Power Measurements 

The amount of electrical power that the heater is using 
is found by placing a voltmeter in series and an ammeter in parallel 
with the power lines returning from the heater. 

The power level is controlled by a field rheostat placed 
in parallel with the DC motor-generator. A 150-ampere circuit breaker 
mmc luded in emne circuit. 

A slip-ring unit must be used to obtain readouts from the 
instrumentation mounted on the rotating heat pipe. A photograph 
of the slip-ring unit is shown in FIGURE 6. It utilizes liquid mercury 
as a viscous contact medium between the revolving and stationary leads. 
This cuts the frictional heating noise to less than 10 microvolts at 
speeds up to 4000 rpm. The slip ring unit has 22 terminals and is 
therefore capable of accomodating 9 thermocouples and 1 pressure 
transducer. The output from the slip-rings is fed to a Hewlett-Packard 
2010 C Data Acquisition System containing an integrating digital 
voltmeter, guarded data amplifier, and guarded cross-bar scanner. This 


system provides printed digital output with an accuracy of +0.5 microvolts. 


Cu PROPOSED OPERATING PROCEDURES 
The operating procedures are divided into three sub-procedures: 
1.) filling procedure, 2.) testing procedure, and 3.) treatment of 


data. 
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1.) Filling Procedure 
Before operating and testing the rotating heat pipe, the 

working fluid must be introduced into the System. This is done in 
such a manner as to keep the amount of non-condensable gases in the 
system to a minimum. These gases reduce the amount of condenser 
surface available to the vapor. Before mete commences, a maximum 
operating saturation pressure is decided upon. A limit on the operating 
pressure is imposed by the amount of pressure the heat pipe is Reiasedu 
to withstand. The maximum operating pressure determines the minimum 
amount of working fluid that is added to the pipe (the maximum is 
determined by the amount of liquid the evaporator annulus is capable 
mG holding) . 

The pipe is placed in a vertical position with the evaporator 
section down and the shaft not in place. The working fluid is added 
and a heat flux is generated by passing a current through the heater 
element. As the liquid evaporates, the vapor drives the non-condensable 
gases to the end of the condenser and out through the hole in the end 
plug. Once it is felt that the gases are eliminated from the system, 
the shaft is screwed into the threaded shaftway in the end plug until 
it seals against a metallic O-ring located on its end. Heat is still 
added and liquid evaporated to determine if there is enough working 
fluid to attain the pre-determined maximum operating saturation pressure. 
Once this is ascertained, the current to the heater is shut off. The 
filling phase of operztion is completed and testing can commence. 

2.) Testing Procedure 

Once filling is completed, the test stand is placed in the 


chosen orientation for the run. Rotation is started and the rpm is 








brought to a pre-determined @leveleior (hewrun es hiss cecheciceoman, 
readings obtained from the optical tachometer and strobe light. 

The spray cooling system is turned on. Finally, a current is passed 
through the heater element. The current is increased in steps with 

a waiting period between each increase to allow the pipe to reach steady- 
state. The readings of the thermocouples, ammeter, voltmeter, quartz 
thermometers, pressure transducer, . coolant pressure tap, and rotameter 
are recorded with each current step. This is continued until the 
maximum saturation pressure level is reached. ‘The same procedure is 

then repeated for a series of increasing rpm values. 

Shut-down is accomplished by first shutting off the current to 
the heater, then decreasing the rpm by use of the remote control unit 
(simply turning the motor off will engage the magnetic disc brake). 
Finally, the coolant system is secured. 

3.) Treatment of Data 

The amount of heat the heat pipe is transporting during 
steady-state operation is the primary quantity to be determined from 
the experimental data. Because of its importance in fixing the 
operating characteristics of the rotating heat pipe, three methods are 
used to obtain the heat flux. 

The first method is to determine a radial temperature profile 
in the evaporator wall from the readings of the thermocouples imbedded 
faeene Wall. dhe heat flux is then obtained from Kourier’s heat 
transfer equation in cylindrical coordinates. 

The heat flux is also obtained by recording the voltage and 
current levels to the heater. Once the boiler heat losses are 
approximated, either by the analytical approach described in Appendix 


A (using different convection coefficients calculated for different 
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speed levels) or the experimental procedure described in ref. [6] , 
the actual heat flux is obtained by subtracting this loss from the 
gross power input. 

The third method for determining the heat flux is by performing 
an energy balance on the condenser cooling section. The coolant 
flow rate is obtained from the flow rotameters and from a calibrated 
curve for the standard nozzles spraying at a fixed pressure, known 
from the pressure tap reading. The coolant's mean temperature is 
measured before and after spraying by quartz thermometers mounted 
in mixing tubes located above and below the spray cooling evilander, 
With this information and knowing the coolant's density and specific 
heat, a heat flux is obtained from a simple energy balance. 

The vapor pressure is another important quantity that is 
monitored during heat pipe operation. It acts as the independent 
variable in the experimental analysis of the rotating heat pipe. ‘The 
vapor pressure is determined both from the reading of the pressure 
transducer in the vapor and the thermocouple mounted in the vapor 
space. The thermocouple gives an indication of the saturation 
temperature of the vapor, which can be converted to an equivalent 
saturation pressure. 

The rotational speed is the adjustable parameter in the 
experimental study of the rotating heat pipe. As indicated previously, 
its value is obtained by both an optical tachometer and electronic 


strobe light. 
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ITI. ANALYTICAL PROGRAM 


py PRELIMINARY ANALYSIS 

In his thesis on the operation of a rotating wickless heat pipe, 
Ballback J7/} discussed the limitations on the rotating pipe due to 
various fluid d‘namic phenomena. Using experimental correlations 
taken from systems approximating the dynamics of the new heat pipe. 
he was able to estimate the -imits of operation imposed by the mechanisms 
Sisentrainment, sonic velocity and critical nucleate boiling. In 
order to obtain a condenser limit for the conically-shaped condenser 
section Ral’ back performed a Nusselt-type analysis for the condensa- 
tion on a truncated cone. Fxisting solutions for rotating condensation 
i, 9] were not applicable to this geometry since they could not satisfy 
the condition of zero condensate velocity at the truncated end of the 
condenser. Ballback!s basic assumptions for such an anlysisS were 
justified by classical theory = experimental correlations. With these 
assumptions and the coordinate system shown in FIGURE 8, he performed 
force balances in the X-and Y-coordinate directions on an infinitesimal 
fluid element. After integration and simplification of these basic 


momentum relationships, an expression for the condensate velocity was 


obtained. 
2 wes 
Ue = [SY YY [ (R+Xsing- 6 cos@) (sind-cos ax )| (1) 
i 
where U, Setluideveliocity (| ft/hr) 


Dp = fluid density (1bm/ft?) 
Bp = fluid dynamic viscosity, (2 bm/ft-hr) 
= rotational speed (1/hr) 
Y = coordinate measuring distance normal to surface, (ft) 
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= condensate film thickness, (ft) 
minimum wall radius in the condenser section, (ft) 


half-cone angle, (rad) 
coordinate measuring distance along the surface, (ft) 


it oli 


§ 
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X 
sf 


slope of the condensate film 





FIGURE 8 Coordinate System for Condensing Limit 
From continuity,an expression was obtained for the condensate mass 


Flow rate. 
3 
n= 27 aa [ (Rptsing ~8e054)(sind-cose9S)]  Kqyrsing) = t*coss (2) 
Ms 


Ballback then assumed that the film is very this (ie, §<< R,) and that 


of «tang (3) 


This assumption allowed the mass flow rate expression, eq'n(2), to 


be simplified to 


i= Aer [(apxemg)?sing] by 


With a steady-state assumption, the rate of condensation equation 


b 9 
ecame h, d=, : dQ, 


gdxX dX 
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(5) 





where Dion latent heat of vaporization of the working fluid, 


(Btu/lbm). 


tt was obtained by considering the heat transfer by conduction 


to an element of surface area. In determining this expression, 
Ballback considered only the thermal resistance due To the condensate 
film and neglected the thermal resistances in the condenser wall and 
outside cooling mechanism. This simplified expression was substituted 
in equation (5) which. after integration, yielded a closed-form solution 
for the condensate film thickness as a function of X. Ballback defined 
a local heat transfer coefficient as the thermal conductivity of the 
condensate divided by its thickness. He integrated this expression 
over the surface of the condenser to obtain a mean heat transfer 
coefficient. With this mean coefficient a total heat rate for the 
condenser was derived. Ballback then coupled this condenser heat rate 
with a conduction equation across the condenser wall to generate a 
condensing limit curve. 

Ballback's results were obtained by neglecting the thermal 
resistances in the Sense wall and the external condenser cooling 
mechanism, and by making the additional simplifying assumption that 


dé<<tang. (3) 
dx 


It is apparent that the above assumption, eq'n(3), is true only for 
large half-cone angles. Also if Ris zero. the mean heat transfer 
coefficient obtained in Ballback's analysis is identical to coefficient 
obtained in an existing boundary layer analysis of condensation on 

a rotating cone performed by Sparrow and Hartnett [9]. They point out 


that their solution is limited to large half-cone angles. 
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Ballback's theory can be improved if it is modified to take 
into account the thermal resistances in the condenser wall and 
condenser cooling mechanism. This is accomplished by including these 
cerms in the incremental heat conduction equation, which, for a thin 
condensate film and condenser wall, becomes 


dQ, 2 7 (REXsinO)(T_-T, 4) 


(6) 
ax 7k +¢/k + 17h 
t W 


where k, = thermal conductivity of the liquid condensate 
(Btu/hr-ft- F) 


kK = thermal conductivity of the stainless-steel 
condenser wall (Btu/hr-ft- F) 

t = thickness of condenser wall, (ft) 

y= 


condenser cooling mechapism convection 
coefficient, (Btu/hr-ft - F) 


Wem this expression is substituted in the rate of condensation equation 
(5) along with Ballback's simplified expression for the condensate mass 
flow rate, a-first-order. non-linear differential equation for the 
condensate film thickness if generated. A closed-form analytical 
solution for&cannot be obtained from this formulation. A Runge-Kutta- 
Gill numerical integration scheme was therefore used to solve eqins(4), 
(5) and (6) using an IBM - 360 Model 67 digital computer and Ballback's 
heat pipe dimensions. A zero initial condensate film thickness was 
assumed to start the integration scheme since this is the only value 
that gives a zero condensate velocity at the truncated end of the 
condenser. Since this theory is still limited to large half-cone 
angles by the neglect of the condensate film slope in the mass flow 
rate expression, the results obtained from the above numerical 
integration are approximate and are assumed to represent an upper limit 


on the condensing limit for rotating heat pipe operation. 
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The condenser limit heat transfer rates are plotted versus 
rotational speed for a saturation temperature of 212°F and half- 
cone angles of one- and two-degrees in FIGURES 9 and 10. The 
outside convection coefficient is a parameter showing the 
dependence of the heat transfer on the condenser surface cooling 
mechanism. With an infinite convection coefficient and a 3° half- 
cone angle, a comparison between Ballback's predicted condensing 
limit and the condensing limit calculated from the modified theory 
is shown in FIGURE ll. The modified theory is seen to predict a 
slightly lower condenser limit heat transfer rate; apparently, 
neglecting the wall ie aieaee during integration leads to higher 


rates of heat transfer. 


B. REVISED ANALYSIS 
The large half-cone angle restriction on the above modified 
analysis can be removed by including the film thickness slope term 


in the expression for the condensate mass flow rate. 


My < mee S| cRpxsmng- §cosd) (sinp-cosg)] [Strpysm 6-5, 5* costlf 


(2) 


When this expression is differentiated and substituted in the rate of 
condensation equation for steady-state operation eq'n(5), a second- 
order differential equation for the condensate film thickness is 


obtained. 
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FIGURE 9 HEAT TRANSFER RATE vs ROTATIONAL SPEED 
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FIGURE 10 Heat Transfer Rate vs Rotational Speed 
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FIGURE 11 Comparison With Ballback's [7] Condensing Limit 
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This expression is valid for all half-cone angles. It reduces upon 
simplification to existing theoretical expressions for a rotating 
disc oe g=90° and for a rotating drum [10] at $=0°. 

A numerical integration method must be utilized to solve 
equations (6) and (7) simultaneously to obtain a solution for the 
condensate film thickness as a function of X and the condenser limit 
heat transfer rate Q,- To start the integration, initial values for 
the heat rate, film thickness, and film thickness slope must be 
determined. Q, is zero at X=0 since there is no area for heat transfer 
to occur (the trumcated base is assumed to be adiabatic). Leppert 
and Nimmo [10] state that for condensation on finite surfaces normal 
to inertial forces, the starting condensate film thickness is a 


function of the minimum film thickness, which occurs at (or very near) 


the surface's edge. This is shown in FIGURE 12. 





a ee 


FIGURE 12 ee of Coniensate Film Under Taf lnence ai Thee 
Forces 
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The minimum thickness depends on the particular overfall condition 
occurring at the surface's edge.Clifton and Chapman[12] use a 
free-overfall edge condition to analyze natural convection on a 
finite horizontal plate. If density variations are replaced by 
rotational acceleration as a driving force, an expression can be 
obtained for the minimum condensate film thickness using a free- 
overfall edge condition. This analysis is sieiened in Appendix C 


with the result that 





oe ace (8) 

Pp (2M* (Rt L, sing)@ 

_where § =condensate film thickness at the overfall, X=L From 
equation (1) in order for the condensate velocity to be zero at 


the trumcated end of the condenser section, the slope of the condensate 


film must initially equal the tangent of the half-cone angle. 
ds - tand. (9) 
3X 


The method of solution of the differential equations for the 
condensate film thickness and condensing limit heat rate resembles 
that for the solution of a system of non-linear differential equations 
with an asymptotic boundary condition. With the initial values of 
the heat rate and condensate slope known, an initial value of the 
condensate thickness is guessed to start the numerical integration 
scheme. With these starting values equation (7) is integrated, 
determining the condensate film thickness as a function of X. The 
film thickness at X=L, is compared to that obtained from the free- 


overfall condition eq'n(8). The overfall thickness is a function 


TBecause of the coordinate system (see FIGURE 8), this is an approx- 
imation for all half-cone angles but zero. The exact boundary condition 
for these angles is not known because the condensate velocity is not 
zero at X=0 
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of the condensate mass flow rate at X=L which is calculated from 
equation (2). If the two thicknesses compare, a solution is 
obtained. If they do not compare, a new guess for the initial film 
thickness is used as a starting value in the integration scheme. 
This procedure is repeated until a satisfactory result is arrived at. 
A Runge-Kutta-Gill numerical integration scheme was used to 
solve equations (6) and (7) simultaneously on an IBM 360 Mod 67 
digital computer. The initial condensate film thickness was iterated 
until the final thickness obtained from the integration and that 
obtained from the overfall condition agreed to within .0004 inches. 
Very favorable results were obtained from this procedure for 
a half-cone angle of zero degrees. However, as the value of the 
half-cone angle was increased, a problem was encountered which 
centered about the instability of the highly non-linear second-order 
differential equation for the condensate film thickness. As is the 
case with more well-known non-linear differential equations with 
asymptotic boundary conditions [14] the instability of the equation 
was highly sensitive to the starting value for the film thickness. 
Consequently, no numerical results were obtained for half-cone angles 
greater than zero. The condensing limit heat transfer rate is 
plotted versus rotational speed for a zero half-cone angle for different 
values of the outside condenser surface convection coefficient, in 
FIGURES 13, 14, and 15. The saturation pressure and therefore temperature 
is a parameter on these graphs. When the thermal resistances in the 
wall and outside cooling mechanism are neglected, the results compare 


with those of Leppert and Nimmo [10] for a rotating cylinder. 
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TV SUMMARY 


A. EXPERIMENTAL GEOMETRY 

i The design and machining of the rotating heat pipe has been 
accomplished as described earlier. It is designed such that visual 
observations can be made of internal pipe mechanisms such as critical 
nucleate boiling and condensate entrainment. The instrumentation system 
is devised so that the parameters of the experimental study such as 

heat flux, saturation pressure, and rotational speed can be monitored 

in more than one manner so that checks on a reading can be obtained. 

A high degree of precision was used in the machining of the pipe 

and the mounting of the instrumentation to keep the degree of un- 
certainity in the determination of the operating characteristics to a 
minimum. 

ce. The ea pipe and drive assembly are mounted on a special test 
stand so that tests can be taken on the device at different orientations 
in the gravity field. 

Dr Safety is insured during operation by enclosing the rotating 
assembly within stainless-steel shielding. Plexi-glass viewports and 
doors in this shielding still allow visual observation of the operation 


to be accomplished. 


i. THEORETICAL STUDY 

a The preliminary theoretical investigation performed by Ballback 
is limited to large half-cone angles and by the fact that the mean 
heat transfer coefficient was obtained by considering the effect of 
the thermal resistance due to only the condensate film and neglecting 


the effects of the condenser wall and outside cooling mechanism. 
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ey After amending Ballback's theoretical results to include the 
thermal resistance effects of the condenser wall and outside cooling 
mechanism, the condenser limit for a particular operating condition is 
seen to decrease (item less heat can be removed than was previously 
expected). 

Dy. Based on Ballback's basic Nusselt-type analysis of condensation 
fea cruncated cone Under the influence of rotational forces, an 
improved theory was obtained that is good for all half-cone angles and 
includes the effects of thermal resistance in the condenser wall and 
outside cooling mechanism. This formulation was coupled with a free 
overfall condition and produced results fora 0° cone angle that are 

in close agreement with those from existing theory for a rotating 
cylinder. Results were not obtained for non zero half-cone angles 
because of the instability of the governing equation, the dependency of 
pais Pendition being on the accuracy of the initial guess for the film 
thickness at the truncated end of the condenser section. 

Hu The condenser limit of the rotating heat pipe is governed by 

the maximum pressure the pipe is designed to withstand and the amount 
of working fluid it can contain. If this peak operating level is 
designed high enough the performance of the pipe may become limited by 
entrainment or critical nucleate boiling before the condenser limit is 
approached. 

5. The performance of the rotating heat pipe is governed by the method 


by which the condenser section is cooled. 
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V. RECOMMENDATIONS 


A. EXPERIMENTAL ASPECTS 

i The heat pipe assembly should be completed and the heat transfer 
capabilities tested for various operating conditions as described in 
wots Study. 

a, Different condenser geometries should be designed, machined, and 
tested. Some of the variations could include the varying of the half- 
cone angle, the addition of some non-wetting material to the condensation 
surface to promote dropwise rather than film condensation, the machining 
of longitudinal grooves on the condensing surface to decrease the effect 
of the interfacial shear on returning condensate, and the placing of 
fins on the outside condenser surface to increase the effective heat 
transfer area. 

oe Different condenser cooling mechanisms, from air to liquid gases 


should be utilized and studied. 


iB. THEORETICAL ASPECTS 

Hs, An iterative scheme such as that described ref. [14] should be 
formulated and included in the computer program to aid in obtaining 
theoretical results for half-cone angles other than zero degrees. 

2 The effect of the vapor should be included in the theoretical 
analysis by solving the continuity, momentum and energy equations in 
the vapor space and coupling this solution with that of the condensate 


by the common shear stress and temperature at their interface. 
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APPENDIX A 


DETERMINATION OF MAXIMUM RADIAL 
HEAT FLUX AND ASSOCIATED LOSSES 


ae. The heater element is made of Chromel A which has a melting 
point of approximately 2u00°F. It is assumed that its maximum 

O 
surface temperature is 2000 F. Using a cylindrical coordinate system 


oe 
a enh T/R, 


where 


R =thermal resistance of stainless steel 
boiler wall (°F-hr-ft/Btu) 
Ine = thermal resistance of inconel heater 
sheath (~F-hr-ft/Btu) 
R,, sthermal resistance of magnesia heater 
MgO O 
dneulation. | f-hr-ft/ Ben) 


The resistance of the fluid annulus is neglected due to boiling. From 
the known geometry of the evaporator section and the known thermal con- 
ductivities of the various elements 


aeons OP nr-ft 
R= "y 5625 =.01237 Btu 
12 


mele 8525 OR-hr-ft 

R05 "7 .8125=.00127 Btu 
1.8325 OF-nr-ft 

1 Sete 

Rye o™ 1.8325=.01275 ~ Btu 
1.2 
R,=.02638°R-hr-ft 
t pam e 
Btu 


Assuming a maximum operating saturation pressure of 3 atmospheres, the 
temperature difference across the resistance is 


T2000-270=1730°F 
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Then 


Cie C0) wee Btu 
Ll. ~ Oaese = 420,000 hr-ft 


The heating surface is approximately 4 inches long and the boiler inside 
surface area is 
pone (2.5625) (4) (1/144)=.275 ft ¢ 


The heat flux is then 


Q 


FEI 76, cy) 
Ay 


hr-ft 





2 


a: The radial heat losses are calculated in much the same manner. 
Neglecting the effect of the thin coating of Sauereisen cement, the 


thermal resistance is 


1 


ae 
Domine? MeO ins Wh 


on 


where 


R, 7 thermal resistance of insulation, (°F-hr-ft/Btu) 


=outside convection coefficient between, the 
insulation and atmosphere(Btu/hr-°F-ft*) 
A=outside surface area of insulation. 

It is assumed that the heat transfer on the outside of the 
rotating evaporator assembly can be approximated by the flow over 
a stationary cyliner. Then ref. [11] can be used to obtain an 
estimate for the outside heat transfer coefficient. 


m 
Nup=C(Rep) for W=3600 rpm 
Rep= 200,000 Re,=(3600) (27 ) (60) (a) 
C= .0239 Vair 
m= .805 
where d=diameter of insulation, 
Re Omeen re cner to 
ee 0 an =kingmatic viscosity of air 
—~ .0277 F-ft-hr/Btu alt (¢t°/hr) 
2800 ¥ 
R, = 1.98=.00116 °F-hr-ft/Btu 
= OmoG 
-=in 1.980_ Oo 
R eg5e = 2OUnOy helt t/ Btu : 
Me Oe ee ee 
iliac 


a1 





inf l2 


R, "2.00 = 6.37 °F-hr-ft/Btu 
nS. ors 


These were obtained by assuming a 3/4-inch layer of insulation with 


a thermal conductivity of .05 Bede ete 


R,=6.41 °F-hr-ft/Btu 
< a 1890 Btu/hr-ft 
2 


GOR 
Cisse a tee Btu/hr-ft 


h can be re-calculated for different rotational speeds, giving a 
relationship between heat loss and rpm. 
De Using the known evaporator geometry and assuming that 19 helical 
turns of heater can be wrapped on the 4-inch heating surface, the total 
length of heater is 

2 nm (1.90625) (19)/127 19 ft 
The resistivity of the heater element is known to be .0785 ohms/ft. 
ies tixes the total resistance in the circuit 

R=(.0785)(19)=1.4887 ohms 
If the maximum voltage that can be supplied is 220 volts, the current 
i 

T=220/1.4887= 150 amperes 
The gross power input is therefore 

Pal R=32.5 kilowatts 
which corresponds to a heat flux of approximately 350,000 Dey eee 
assuming losses of 10%. 
yy It is therefore seen that the total heat flux capable of being 
conducted across the evaporator wall is limited by the amount of heat 


that can generated by the heater. 
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APPENDIX B 
PRELIMINARY UNCERTAINTY ANALYSIS 


Assume that the heat flux is to be determined at a surface at 
a radial distance R by the readings of two thermocouples placed at 


e e ° 4 e ® e 
radial positions ro and ree The governing equation is 


an Kies 
VT Ap R In(r/t7) 
-where 
T= temperature reading obtained at i. 
T, = temperature reading obtained at ra 
A, = surface area at the radial position R 
k = thermal conductivity of material being used 
- dq= 28 4 K+29 Edt +59. oF, br, +33. dr, 2438, dr? + S4dr 
“From the governing equation 
i PRG) +7 Rhel 
ye dq _ “KET, ) 
ae RR * Bin a/e) 
-K(CT. -T; ) Kea 


Q 
ie, “WR DoGee/e] > aR TERT 


. 1t is seen 


a dT. oT; off Or, 
i een) R r’ ri, In(r,/n’) i r Te i) 
: it is assus:d that all the uncertainties are additive and that the 


values of k and R are known with high accuracy, it is seen that the 


uncertainty in the determination of the heat flux depends on the. 
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uncertainty in the readings of the thermocouples and in their positioning. 
It is also seen that the farther apart the thermocouples are placed, 
the larger the value of the temperature difference and the logarithmic 


terms become and the smaller the effect of their associated uncertainties 


become. 
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PereND Dae 
DERIVATION OF PRE 
OVERFALL END CONDITION 
From ref. fi3}it is seen that during flow down a mild, inclined slope, 


critical flow will occur at a sudden change in slope as shown below. 





From energy considerations of the system 


ay (4 pr? 


ax (1-Fr') = S. ~ Sp 
where 
Fr = Froude Number a 
/bS 
S._ = bed slo meee 
oun eee Ox 
y* 
5 =< ricrtren slope = —— 
f 2 
c R 
Ss 
R. = wetted mean radius (=A/P) 
ec = Chezy constant (similar to friction factor, f) 


In the region of mild slope, SPs. Since the slope is so mild. On the 


other hand, in the region of steep slope, S <5): Therefore there 


f 


must exist a point, or region, where Sp = So" neice 


2 


dy = 
a (1-Fr') = 0 
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wince the slope of the free surface 1s ne vs2cre eae lew a ela yee 


Vis pOsne 
Fr = 1,0 


which is the critical condition. Thus, 


V/ vey, = 1.0 
or 
2 
Vere 


where y. is the crivtvreaiadep vit 


For flow over a sharp corner, the flow can separate from the 


steep wall as shown below. 





a 
7 
YA 
ie 
Va 


This is called a free overfall. In this case, there is experimental 


data [15] to show that 


Y 
sae oP Lye 


The value of Ue occurs very near the edge of the plate, so that in the 


case of the condensate film flowing on the condenser wall 


O.7 y at x = L, 
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where 


§ 


a - condensate film thickness at the overfall 


Le = the condenser length 
2 
Yoo Weare 
Here, V is the average velocity of the condensate given by 


Va } {3 a 4 te 
a u TY Y SC ——— 
AMR § 0 # ZiTR 6 Pg 
where 


Up = velocity of the condensate 


condensate mass flow rate 


ie 
Per 


Also, the acceleration of gravity, g, is replaced by the rotational 


density of the condensate 


acceleration term 
2 
g= 0 RK 
3 a (R +X sin) 
where : 


rotational speed 
small radius of the condenser cone 
half-cone angle 


QW E 


Therefore, 2 
™m 
¥y = 


y gs w(Ri+X sing) (anR§ )* 


or, at X=L 
Cc e2 
0.7 Me 


0  ——E——E— 
be phew X(7n)?(R +L sin)? 


Finally, as the free overfall condition in the condenser, 





0.7 yn 
C1) ptw? (zm)? (R +L sing) 





ey 
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